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Abstract 

We have theoretically investigated the properties of electronic transport in graphene heterostruc- 
tures, which are consisted of two different graphene superlattices with one-dimensional periodic 
potentials. It is found that such heterostructures possess an unusual tunneling state occurring 
inside the original forbidden gaps, and the electronic conductance is greatly enhanced and Fano 
factor is strongly suppressed near the energy of the tunneling state. Finally we present the match- 
ing condition of the impedance of the pseudospin wave for occuring the tunneling state by using 
the Bloch-wave expansion method. 
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I. INTRODUCTION 



Graphene has attracted a lot of attention since its discovery in 2004 [1]. The linear 
dispersion property near the Dirac point leads the low-energy charge carriers obeying the 
massless relativistic Dirac equation. Thus many unique electronic and transport properties 
such as half-integer quantum Hall effect and Klein paradox (4J are demonstrated in 

graphene. 

Currently, most investigations are focused on the properties of graphene itself. There 



are only a few works on graphene-based heterostructures, 
heterostructures 



or example, gated graphene 

Baphene nanoribbon heterostructures [6], and graphene/hexagonal- 
]. Heterostructures are usually referred to the combination of two 
different materials/structures [8|. According to the analog between the one-dimensional 



photonic crystals containing left-handed metamaterials |9Hllj and the graphene-based one- 



dimensional periodic superlattices 



many optical properties may be introduced into 



the graphene-based electronic systems. In this paper, we will theoretically study the 



graphene heterostructures consisted of two different graphene-based super 



attices (GSs), 



which can be constructed by electrostatic |l5H20j or magnetic barriers [2l|, |22|. Since the 
electronic pseudospin wave may be highly localized at the interface between two different 
GSs, it would be expected that such heterostructures may possess an unusual tunneling 
mode inside the zero-averaged wave- number gap [k| . The location of such a tunneling state 
satisfies the impedance matching condition of the electronic pseudospin wave. 



II. THEORETICAL FORMULA 



Now let us consider an electron of energy E and wavevector ko = EjUvp (with the 
Fermi velocity vf ~ 10 6 m/s) be incident at an angle 8q onto a graphene heterostructure 
(AB) n (CD) m , as shown in Fig. 1, where A, B, C and D indicate different potentials, and n 
and m are the period numbers. As we know that the Hamiltonian of a low-energy electron 
in graphene is given by Hj = vfo ■ p + Vj(x)I, here p = (—ihJ^, —ih-^) is a momentum 
operator, a = (a x ,a y ), a x and a y are Pauli matrices, I is a 2 x 2 unit matrix, and Vj(x) 
is the jth electrostatic potential. This Hamiltonian acts on a two-component pseudospin 
state ^ = (ipi,ip2) T '> where ip x and ip 2 are t wo components of the pseudospin wave (here 



2 



we omit the term e tkyV due to the translation invariance in the y direction). Therefore, the 
wavefunctions ip l2 at x and x + Ax inside the jth potential can be related via a transfer 
matrix: \\\ 

cos(k x jAx—9j) ■ sin(fcj; jAx) 

Mj(A*. E, k,) = | "'^ | , (1) 



where k x j =sign(kj)*/kj — k 2 , is the x component of the wavevector kj = [E — Vj(x)]/hvF 



for kj > ky, otherwise k x j = iyk 2 — k 2 . Using the boundary conditions, both the reflection 



and transmission coefficients could be obtained by the transfer matrix method, 



13| 



r(E,ky) = 7 / — v 1 , (2) 

. fJ? , \ 2cosg 

where so = e l9 °, s e = e l9e , 6 e is an exit angle (see Fig. 1), and Xij are the el- 

(\ 2n+2m 
Xn X\2 \ -p-i- 
= II Mj{Wj, E,ky) from the 
3^21 3^22 I j=l 

incident to exit ends. Once the coefficients r and t are obtained, the total con- 
ductance G of the system at zero temperature can also be calculated by using the 
Biittiker formula 23], G = Gq f^ 2 T(E, k y ) cos 8od9 , where T = \t(E, k y )\ 2 is the 
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transmittivity and Gq = 2e 2 mvFL y /h 2 , and L y is the width of the graphene strip in 



the y direction. Furthermore, the Fano factor in this system is given by |24( F = 
T(l — T) cos9od9oj / ( f^ 2 T cos 6 od9o\ Therefore the reflection, transmission, 
conductance, and Fano factor inside such graphene-based heterostructures can be obtained 
by numerical calculations. 



III. NUMERICAL RESULTS 



Based on the above formulae, we can calculate the changes of the transmission, conduc- 
tance and Fano factor of an electron passing through the heterostructures. In Fig. 2 and 
Fig. 3, we show the typical electronic transmission, conductance and Fano factor for different 
situations with the parameters n = 9, m = 16, Va = lOmeV, Vb = 80meV, Vc = HOmeV, 
Vd = — 20meV, and wa = wb = wc = wd = 20nm. For a single GS (AB) n or (CD) m , from 
Fig. 2(a) and Fig. 3, the forbidden gaps do exist near the new Dirac point at E — 45meV, 
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which is called the zero-averaged wave- number gap in Ref. [13||, and the corresponding con- 
ductance is minimal and the Fano factor equals to 1/3 at the new Dirac point of E = 45meV 
for both cases of (AB) n and (CD) m . Here we would like to point out that the condition 



for appearing the Dirac point inside the GSs are discussed in detail in Refs. |12Ml4J|. For a 
special case with wa = »b and wq = wd, the condition is simply given by E = (Va + Vs)/2 
and E = (V C + Vb)/2, derived from Eq. (41) in Ref. (l3(] . 

However, for the heterostructure (AB) n (CD) m , see Fig. 2(b), a new tunneling mode 
occurs inside the forbidden gap near the Dirac point. Due to the existence of this tunneling 
mode, compared with the case of the individual (AB) n or (CD) m , for the heterostructure 
the electronic conductance is greatly enhanced and the Fano factor is strongly suppressed, 
as shown in Figs. 3(a) and 3(b). In fact, this new tunneling mode is due to the formation 
of the highly localized state at the interface between (AB) n and (CD) m , which could be 
clearly seen from the density distributions of the electronic pseudospin wavefunctions in Fig. 
4. The probability densities are highly localized at the interface between (AB) n and (CD) m , 
which is denoted by the vertical dot line. 

IV. MATCHING CONDITION OF THE TUNNELING MODE IN FORBIDDEN 
GAPS 

In order to obtain the condition of this new tunneling mode, we use the so-called matching 
condition of the impedance of the pseudospin wave. In infinite periodic barriers, y obeys 
the Bloch-Floquet theorem. By using the Bloch-wave expansion method in Ref. [25|, inside 
the GS (AB) n or (CD) m , * A b,cd can be written as 

=H*y ^ j e^ x + H T L ^ ^ j e~^ x , (4) 

where (i p x is the x component of the pseudospin Bloch-wave vector and satisfies the eigen 
equation: Q p ^f p = e ±J/3 * Ap \l/p [Q p is the matrix of the unit cell for (AB) n or (CD) m with 
p = AB or CD], and A p is the length of unit cell. (1,Z±) T are the eigenvectors of Q p 
corresponding to the eigenvalues e ±l ^ Dp J and H± are the amplitudes of the forward (+) 
and backward (-) Bloch waves. From Ref. [13J, \& at both the incident and exit regions 
can be, respectively, expressed by ^(CT) = + r) (1,£ _) , and ^(Xg) = t$j (1,^+) , 
where £ - = (e t9 ° — re~ ie °)/(l + r) and ^ a ,+ = e lde are the impedances at the incident and 



exit ends, and $j is the incident wavefunction. Inside the GS (AB) n } from Eq. fl4]), we can 
define \I/ab(0 + ) = H Q + (1,£ +) T and $ ab{x~d) = H x - fl)^-) , where £ + is the impedance 
at x = + and f^.- is the impedance at the interface x = x^ (xd = tiA^b), and H +(H X -) 
is a coefficient. Therefore, we can obtain 

yABfc yAB\\2 yAB(c yAB\ 

p _ z + Uo+ ~ z - ) A AB — L ~ l£0+ ~ Z + J /r\ 

^ (e + - z AB )\\ B - (^ 0+ - z?) 5 u 

with Aab = e m ^ BKAB . Similarly, the impedance of the pseudospin wave at x = inside 
the GS (CD) m is given by 

p _ z c + D (^-z^)-z^(^-z^)M D 

where Xcd = e tm,3x Ac -°. For a perfect tunneling mode in heterostructures, when and only 
when the impedances of the pseudospin waves at two sides of the interface between (AB) n 
and (CD) m satisfy 

ic--e 4 i = o, (7) 

the zero reflection can be exactly reached. In this case, the reflection coefficient is also zero, 
i. e., r = 0, thus £ + = £ - = e ld ° at the incident end, and £ x - = £ x + = e l9e at the exit 
end. For the imperfect tunneling process, the condition of Eq. (j7j) should be replaced by 
lCc~ — £x+\ ~^ m i n - m Fig. 2(b), we plot the curve for the value of \£ x - — £ x +\, see the 
dashed line. From Fig. 2(b), it is clear that when the value of \£ x - — ^+| is close to zero 
(within the original forbidden gaps), there is a tunneling peak with near- unity transmittance. 
According to the analog between one-dimensional photonic crystal heterostructures and 
the graphene-based heterostructure consisted of two different superlattices, we could safely 
say that in our case the tunneling state is the electronic Tamm state. It should be mentioned 
that, for the other zeroes of \£ x - — £ x +\ (which are located outside of the original forbidden 
gaps), they correspond to the resonances of the propagation states, rather than the localized 
states. 

So far we only consider the GSs consisted of square barriers and wells. In fact, the 
above results about the tunneling state will also occur inside the graphene heterostructures, 
which are consisted of other periodic potentials such as cosine or sine potentials [3]. For 
a comparison, in Fig. 5, we also plot the electronic transmittance for the heterostructures 
consisted of different sine potentials. It is clearly seen that the new tunneling state also 
appears inside the original forbidden gaps under the optimal period numbers. 



V. CONCLUSION 



In summary, we have investigated the electronic transport inside the graphene-based 
heterostructures consisted of two different GSs. It is shown that in such heterostructures 
there is an unusual tunneling state inside the zero-averaged wave-number gaps, which leads 
to the enhancement of the electronic conductance and the suppression of the Fano factor. 
The reason of occurring such a tunneling state is due to the highly localization of the 
electronic pseudospin wave at the interface between two GSs. Finally the matching condition 
of the impedance of the pseudospin wave is analytically presented. The parameters used 



in our study are technically feasible as demonstrated in experimental studies 26|-|29). For 



example, in Ref. 26|, superlattice patterns with periodicity as small as 5 nm have been 
imprinted on graphene using the electron-beam induced deposition. The phenomenon for 
occurring such a tunneling state in graphene heterostructures is potentially useful in filtering 
single-e nerg y electron and is also hopefully to use in the design of graphene-based electronic 



devices 



30]. 
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Figure Captions 



Fig. 1. (Color online) (a) Schematic of a heterostructure consisted of two different 
graphene superlattices, (AB) n and (CD) m . (b) The schematic profiles of the potentials. 

Fig. 2. (Color online) Electronic transmittances through (a) the individual graphene su- 
perlattice (AB) n or (CD) m and (b) the heterostructure (AB) n (CD) m , with 6 = 10°. In 
(b), the dashed line denotes the value of |£ x - — £ x +|- The parameters are V A = lOmeV, 
V B = 80meV, V c = HOmeV, V D = — 20meV, wa = w B = w c = w D = 20nm, n — 9, and 
m = 16. 

Fig. 3. (Color online) Changes of (a) conductance and (b) Fano factor in different 
structures. Other parameters are the same as in Fig. 2. 

Fig. 4. (Color online) Changes of probability densities |V , 12 | 2 inside the heterostructure, 
at the energy of the tunneling mode, E rs 44.688 meV, and #o — 10°. Other parameters are 
the same as in Fig. 2. 

Fig. 5. (Color online) Electronic transmittances through different structures with the 
incident angle 6 = 10°: Case 1, V AB (x) = 45 - (35vr/2) sin(7ra/20)meV for < x < 360nm; 
Case 2, V CD {x) = V CD (x) = 45 + (65vr/2) sin(7rx/20)meV for < x < 1080nm; Case 3, the 
combination of Cases 1 and 2. 
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FIG. 1: (Color Online) (a) Schematic of a heterostructure consisted of two different graphene 
superlattices, (AB) n and (CD) m . (b) The schematic profiles of the potentials. 
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FIG. 2: (Color online) Electronic transmittances through (a) the individual graphene superlattice 
(AB) n or {CD) m and (b) the heterostructure (AB) n (CD) m , with O = 10°. In (b), the dashed line 
denotes the value of |£ - — The parameters are Va = lOmeV, Vb = 80meV, Vc = HOmeV, 

X D D 

Vd = — 20meV, wa = wb = wc = wd = 20nm, n = 9, and m = 16. 
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FIG. 3: (Color online) Color online) Changes of (a) conductance and (b) Fano factor in different 
structures. The structral parameters are the same as in Fig. 2. 
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FIG. 4: (Color online) Changes of probability densities \ipi 2 | 2 of the pseudospin wave inside the 
graphene-based heterostructure, at the energy of the tunneling mode, E ~ 44.688 meV and 9q = 
10°. The vertical dot line denotes the interface between (AB) 9 and (CD) 16 . The other parameters 
are the same as in Fig. 2. 




FIG. 5: (Color online) Electronic transmittances through different structures with the incident 
angle O = 10°: Case 1, V A b(x) = 45 - (35vr/2) sin(vra;/20)meV for < x < 360nm; Case 2, 
Vcd(x) = Vcd{x) = 45 + (657r/2) sin(7rx/20)meV for < x < 1080nm; Case 3, the combination of 
Cases 1 and 2. 
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